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Abstract We discuss how quantum dimer models may be used to provide "proofs 
of principle" for the existence of exotic magnetic phases in quantum spin systems. 
The material presented here is an overview of some of the results of Refs. lH and 

ii. 



1 Introduction 

Consider a system of quantum spins on a lattice with antiferromagnetic interactions. 
The energy of a given pair of spins is minimized by forming a singlet but as a spin 
can form a singlet with only one other spin, the system is unable to simultaneously 
optimize all of its interactions. In this sense, the system is said to be frustrated. 
The most common solution to this problem is for the spins not to form singlets but 
instead to magnetically order, as in the classical Neel state. However, one may also 
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envision non-magnetic states where the order is contained in how the singlet pairs 
organize. In this paper, we concentrate on phases where the singlets form between 
nearest-neighbor spins. Such phases are relevant in systems having a spin gap. 

A nearest-neighbor singlet, also called a short-range valence bond, may be repre- 
sented by drawing a dimer across the corresponding link of the lattice. Arrangements 
of singlet pairs correspond to dimer coverings of the lattice, where the hard-core 
condition of one dimer per site captures the frustration. This picture may be ab- 
stracted to a low energy description where the fundamental degrees of freedom are 
the dimers themselves. Quantum dimer models IT] [2]| describe how exotic phase di- 
agrams can arise from the competition between quantum fluctuations of the dimers 
and various potential energies and are, perhaps, the simplest models containing the 
physics of frustration. 

The suggestion that microscopic frustration may lead to nontrivial phases such 
as incommensurate crystals |[3] or liquids of resonating singlets j4] is a rather old 
one. The topic is of current interest due to proposals that such phases may be 
relevant to understanding high Tc superconductors and other strongly correlated 
systems Is] mi?]. In this context, we may ask how a particular exotic phase can, in 
principle, arise from a microscopic Hamiltonian that is local and does not break any 
lattice symmetries. One way to answer such questions of principle is to construct an 
explicit model having the phase in question. Quantum dimer models are well suited 
for this purpose because of the relative simplicity of the dimer Hilbert space. This 
leads to the question of how these effective models can arise from models with more 
physical degrees of freedom, such as SU(2) invariant spin models. 

In section m we review some facts about the simplest dimer models. In section 
[3] we outline the construction of a dimer model that contains a devil's staircase 
of crystalline phases, including states with arbitrarily long (and hence incommen- 
surate) periods. In section |4] we discuss how dimer models may be mapped onto 
SU(2) invariant spin models on decorated lattices. In particular, the existence of a 
liquid phase of resonating dimers in the simplest triangular lattice dimer model im- 
plies the existence of a liquid phase of resonating singlets on a decorated triangular 
lattice. This paper gives an overview of some of the results in Refs. ||8l and ||9l. 



2 Basic models 

The Hilbert space of quantum dimer models is defined by taking each dimer cov- 
ering of the lattice as a basis vector The inner product is usually defined by taking 
different dimer coverings to be orthogonal. The dimer Hilbert space may be par- 
titioned into topological sectors labeled by a pair of winding numbers (Figure [U. 
The winding numbers are global invariants in that they are not affected by local re- 
arrangements of dimers. In particular, local operators (including the Hamiltonian) 
will not have matrix elements between states in different sectors. For bipartite lat- 
tices, the number of sectors is extensive while for non-bipartite lattices, there are 
four sectors corresponding to each winding number being either even or odd. 
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Fig. 1 Winding numbers for the square lattice case. The thick horizontal line extends around the 
lattice. If the vertical lines of lattice are labeled alternately A and B, then the winding number is 
defined as N/^ — Ng where N/^ g is the number of dimers crossing the thick line along A,B lines. For 
the section of the lattice shown, N/^ = 3 and Ng = 0. One may verify that local dimer rearrangements 
will not change this number. In the triangular lattice case, the relevant quantity is defined by just 
counting the number of dimers intersecting this reference line. Whether this number is odd or even 
is a topological invariant but the value of the number can change with local rearrangements of the 
dimers. 



The original quantum dimer model of Rokhsar and KivelsonQ was defined on 
the square lattice with the Hamiltonian: 

H = 2-t(| I iXn I +h.c.)+v(| I iXi 1 1 +1 nX n I ) 

QDM 

The two terms are projection operators that flip (the t term) or count (the v term) 
the dimers if the plaquette in question has parallel dimers, which we call a "flip- 
pable" plaquette, and annihilate the state otherwise. The sum is over all plaquettes 
in the lattice. The model may be generalized to other lattices where the flippable 
plaquettes now correspond to the dimers occupying alternate bonds of the minimal 
(even-length) resonance loop in the lattice. For the triangular lattice ifTOl, these loops 
are still length 4 but there are now six distinct ways of having parallel dimers instead 
of just two for the square lattice. For the honeycomb lattice, the loops are length six 
and the projection operators are three-dimer moves analogous to the benzene res- 
onance. For the pentagonal lattice|81, the loops are length eight, the operators are 
four-dimer moves and so on. 

Models of the Rokhsar-Kivelson type have certain generic properties. When 
V > f, Eq.[I]is positive definite so the system seeks to minimize the number of flip- 
pable plaquettes. In particular, many lattices, including the ones mentioned above, 
may be covered without having any flippable plaquettes and such "staggered" states 
will be zero energy ground states in this limit. When — v ^ f , the system seeks to 
maximize the number of flippable plaquettes so selects an analog of the "columnar" 
state. When v = f , Eq. [T]is again positive definite, but in addition to the staggered 
states, there are a number of liquid states which also have zero energy [2 1. The liq- 
uid wavefunctions may be written as equal amplitude superpositions = Y,c k) 
where the sum is over all states that may be connected by repeatedly applying the 
flip (kinetic) term in Eq.[T] There will be at least one such state in every topological 
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sector. For the square lattice, there is only one state in each sector because the flip 
term is believed to be ergodic within the sectors but additional subtleties may appear 
in other lattices jTOl . Because these are zero energy states, any combination of them 
will be equally valid as a ground state, including the equal amplitude superposition 
of all dimer coverings. 

The liquid states at v = f, which is called the RK point, have no local order 
but they do carry the global quantum numbers of their topological sectors. How- 
ever, the nature of the liquid states depends intimately on the lattice geometry. For 
some bipartite lattices, including the square and honeycomb, these states have al- 
gebraically decaying correlations, are gapless, and the RK point is a critical point 
separating crystalline phases with minimal and maximal winding number. Field the- 
oretic studies ifTTl fT2l suggest that perturbations about this critical point can stabi- 
lize crystalline phases with intermediate winding number and, generically, a devil's 
staircase of such states will occur in between the states of minimal and maximal 
winding. In particular, it was noted that incommensurate crvstals lfTSl could be ob- 
tained by suitable tuning of the perturbation. 

For some non-bipartite lattices, including the triangular and kagome, the correla- 
tions decay exponentially, the liquid states are gapped, and the RK point is the edge 
of aliquidphasef lOl. This was the first example of a model containing a liquid phase 
of resonating bonds and validated earlier suggestions of a physical mechanism by 
which such a phase could be stabilized, namely geometric frustration enhanced by 
strong ring exchange| 14]. For the kagome lattice, a simpler dimer liquid phase, with 
exactly zero correlation beyond one lattice spacing, was obtained in Ref. ifTsl using 
a different construction with only flip terms. 

In the next section, we discuss how a different kind of exotic structure can be 
realized in a dimer model. Then, in section|4] we show how these proofs of principle 
in dimer models may be extended to SU(2)-invariant spin systems. 



3 Modulated states and the devil's staircase 

We now outline a dimer model that shows a devil's staircase of crystalline phases 
where the period becomes arbitrarily large |9|. The following construction is valid 
in the limit of strong coupling and weak (quantum) fluctuations, in contrast to the 
RK point which corresponds to the opposite limit of strong fluctuations. In this 
sense, the connection between the following devil's staircase and the one predicted 
by field theory to occur near bipartite RK points (mentioned above) is tempting but 
speculative. Similarly, we work on the square lattice for convenience but we do not 
believe this choice is so important in the strong coupling regime. 

Our strategy is to perturb a model that contains a degenerate point separating 
crystalline states of different winding number In particular, we consider the follow- 
ing diagonal Hamiltonian: 
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i-a(iexBI -bpxni + ciegxggK dlfflxffll 

+ ....) 

(2) 

where the dots denote terms related by lattice symmetries to ones shown, which 
will appear with the same coefficient. Note that Eq.|2]is local and does not break lat- 
tice symmetries. The coefficients satisfy p,q > c,d > a,b > but fine tuning is not 
required. Terms a and b are competing attractive interactions while the remaining 
terms are repulsive. This Hamiltonian is designed to favor states where every dimer 
participates in only two attractive bonds of the same type. Terms c and d penalize ar- 
rangements with dimers participating in more than two attractive bonds |fT6l . Terms 
p and q serve a technical purpose discussed in Ref. ||9l- 

The quantum system described by the Hamiltonian in Eq.|2]has the zero temper- 
ature phase diagram shown in Fig. [2] As the figure shows, the system prefers exclu- 
sively a or: b bonds depending on which coefficient is larger. When a ^ b, there is 
a large degeneracy where the preferred states involve thin staggered domains sepa- 
rating columnar regions of opposite orientation; the staggered domains come in two 
orientations as illustrated in the figure. 



+ P 



+ q| 



X 




columnar ; herringbone 



(highly degenerate) 




Fig. 2 Ground state phase diagram of the parent Hamiltonian Ho as a function of the parameter 
a — b (see Eq.O- When a — b — 0, there are a number of degenerate ground states. The maximally 
staggered configuration is commonly called the "herringbone" state. 



We perturb this model with a resonance term that is equivalent to two actions of 
the flip term in Eq. [T] 

-tv = 2-t(ig xBl+lcnXm l+h.c.) 



6 



K. S. Raman, E. Fradkin, R. Moessner, S. Papanikolaou and S. L. Sondhi 



where < f <C fl,^. We now describe how the existence of the devil's staircase fol- 
lows by perturbing in t at successively higher orders. Before doing this, we remark 
that the main reason for choosing a somewhat complicated form for Eq.|2]and a non- 
standard resonance term (i.e. Eq.O is that these choices simplify the calculational 
details. As discussed in Ref. ||9l, we anticipate that this construction can be made to 
work for simpler diagonal Hamiltonians and with the familiar two-dimer resonance. 

Referring to Fig.lH we see that the operator[3]does not affect the phases on either 
side of the degenerate point but will partially lift the degeneracy ata^b. In partic- 
ular, this term causes the staggered domain walls to fluctuate which stabilizes them 
at second order in perturbation theory. In this sense, our model is a quantum analog 
of the Pokrovsky-Talapov model of fluctuating domain walls 1 171 . In that classi- 
cal model, thermal fluctuations give entropy to the domain walls and the competi- 
tion between this and their energy cost induces a commensurate-incommensurate 
transition! 3 il to a striped phase where the walls arrange periodically with spacing 
that can be large compared to the interaction scale. In the present case, quantum 
fluctuations play the role of temperature and cause modulated phases to appear. 

A technical issue is that the extent to which a staggered domain is stabilized 
by second order processes depends on its environment. One consequence of this is 
the favoring of states, such as the [In] sequence (Figure [3]i, where the staggered 
domains have the same orientation. Also, while the [11] state clearly has the max- 
imum number of resonances, these resonances are individually weaker than those 
in the [12] state because of the higher energy virtual states that are involved (due 
to term c in Eq. |2]i. If the [11] phase is selected when a = b, then by slightly de- 
tuning from this point (by an amount of order t^), we will reach a value of b where 
the "self-energy" of a columnar line jTSl becomes degenerate with the energy of 
a resonance-stabilized staggered line in the [11] state. However, b must be further 
increased before the advantage of [12] staggered lines is similarly compensated. 
Therefore, both the [11] and [12] phases appear in regions of width ^ between 
the initial phases and to second order in r, we obtain the phase diagram in Fig. [5^. 
Note that the winding number increases as we move from left to right in the phase 
diagram. 




(a) (b) (c) 



Fig. 3 The [11], [12], and [13] states involve staggered domains separated by one, two, and three 
columnar units respectively. In these states, the staggered domains have the same orientation. We 
can similarly define [14], [15] etc. 

On the [1 1]-[12] and [12]-columnar boundaries, the phases and states with inter- 
mediate winding are degenerate to order f ^ but these degeneracies will be lifted at 
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higher orders. Considering the [12]-columnar boundary at fourth order in perturba- 
tion theory, the relevant resonances are shown in Figure |4] Process (a) destabiHzes 
the [12] state due to the high energy virtual state involving term q in Eq. [2] (there 
is an analogous process involving term p) and its contribution is larger than pro- 
cess (b) which is stabilizing. In the [13] state, the analog of process (a) will not 
contribute because the flipped clusters are disconnected which means that the en- 
ergy and wavefunction terms in the perturbation series exactly cancel. However, 
process (b) will occur and stabilize a [13] phase in a region of width ^ between 
the [12] and columnar states. The argument may be applied inductively. Along the 
[In] -columnar boundary, at (2n)th order in the perturbation, the analog of the Figure 
will destabilize the [In] state relative to states of lesser winding while the analog 
of Figure lib selects the [l,nH-l] state from this set. The new [l,nH-l] phase will oc- 
cupy a region of width ~ t^" in the phase diagram and so on. While higher orders in 
the perturbation theory involve increasingly complicated resonances, including fluc- 
tuations causing the staggered lines to effectively "break apart", the calculation is 
designed so that these terms amount to self-energy corrections that simply move the 
boundaries. The stabilization of phases in the [In] sequence will always be governed 
by the analogs of the straight resonances in Figure|4] 



initial state 




Fig. 4 The fourth order resonances driving the transition. We emphasize that these figures rep- 
resent terms that appear in fourth order perturbation theory, not additional terms added to the 
Hamiltonian. The circled cluster in (a) refers to term q in Eq.|2] 



The [11]-[12] boundary can similarly open at higher orders and in Ref. ||9l, we 
verify that a [11-12] phase (the notation means that the repeating unit is a stag- 
gered domain followed by one columnar unit followed by another staggered do- 
main followed by two columnar units) is stabilized at sixth order between the [11] 
and [12] states. The finer boundaries can also open and we have also verified that a 
[1 1 — (12)^] phase occurs at eighth order between the [11-12] and [12] phases. We 
are less certain about this fine structure because the resonances involved are more 
complicated than in the primary [In] sequence. However, we may speculate that the 
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opening of boundaries will continue to finer scales, at least for some range of pa- 
rameters. In this sense, the phase diagram will most generically be described by an 
incomplete devil's staircase, as shown in FigurejSj)- Moving from left to right in Fig- 
ure 5, the system passes through a series of crystalline states with increasing wind- 
ing number via first-order phase transitions. We note that this phase diagram and the 
structure of the calculation is similar to the 3D classical ANNNI model lfT9l l20l 1211 . 




Fig. 5 (a) Ground state phase diagram of H = Ho + tV from second order perturbation theory. The 
width of the [11] and [12] phases are order t^. (b) Sketch of the ground state phase diagram to all 
orders in perturbation theory. While the opening of the [1 1]-[12] boundary is explicitly indicated, 
the other boundaries will also open. The collection of phases forms a devil's staircase. Note: These 
figures are not the result of a simulation but a sketch of the general properties of the phase diagram. 



4 SU(2)-invariaiit realizations 

Having discussed some examples of how exotic phases can arise in dimer models, 
we now turn to the question of what this implies for systems with more realistic 
degrees of freedom. The interpretation of a dimer as a nearest-neighbor singlet or 
valence bond involves two simplifications. The dimers have no orientation and dif- 
ferent dimer coverings are orthogonal by definition. In contrast, valence bonds are 
oriented II22II and different valence bond coverings are not orthogonal. There is also 
the issue of projecting the spin Hilbert space onto the much smaller dimer Hilbert 
space. 

In this section, we discuss one way in which these issues may be resolved allow- 
ing the dimer model phase diagrams to be realized in SU(2)-invariant systems. We 
concentrate on the square lattice but the arguments may be generalized to other lat- 
tices and higher dimensions. The interested reader may consult Ref. lH for details. 
We would also like to point out a complementary approach by Fuiimoto ||23l . 

We construct the spin model on a modified square lattice where the links are dec- 
orated with an even number of additional sites (Fig.|6l). On each lattice site ;, we 
define an operator hj that projects the cluster of that spin and its (two or four) neigh- 
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bors onto its highest spin state ll24l|25ll26l (either 3/2 or 5/2). The parent Hamilto- 
nian is a sum of these operators, one for each lattice site: 

Ho = Y,Oiihi (4) 

/■ 

where a, is a positive constant that, in principle, may vary with i. The SU(2)- 
invariance of Eq. |4] may be seen by explicitly writing down the operators. Re- 
ferring to the figure, = 5^ — {\){\) where S = S^i +Si +Se and hkewise 
for the other link sites. Similarly, In = [S^ - {\){l)][S^ - (|)(§)] where S = 
Si + Sal + Sfoi + Sti + S^/1 and likewise for the other corner sites. 



1 1 





'c1 






d1 \ 


' a1 


a2 " 






!l ^ 


2 


.b2 


b3' 


3 


4 


.b4 



a3 a4 



Fig. 6 Decorated square 
lattice for the case N = 4. 



As these expressions indicate, /i, is a product of terms that annihilate the lower 
spin sectors so only the highest spin state remains. Among the configurations anni- 
hilated by hi is the case where two spins in cluster / form a singlet. This is because 
the highest spin state must be symmetric under interchange. A consequence of all 
this is that any configuration where every spin is in a singlet with one of its neigh- 
bors, which we call a valence bond state, will be annihilated by Eq.|4] Because Hq is 
a sum of positive definite operators, its eigenvalues are positive so such annihilated 
states will be zero energy ground states. Even a minimal decoration of N = 2 en- 
sures that the valence bond states (and their superpositions) are the only zero energy 
states. 

In the valence bond states, the chains forming the links of the decorated square 
lattice are in one of two possible dimerizations. One of the dimerizations involves 
the sites of the original square lattice while the other involves only the decorated 
sites. In the original square lattice, we may represent the former case by drawing a 
dimer across the link and the latter by an empty link. Therefore, the valence bond 
states correspond exactly with dimer coverings of the lattice but an important tech- 
nical difference is the non-orthogonality of the valence bond basis. However, the 
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decoration ensures that the magnitude of the overlap between different valence bond 
states will always be exponentially small in the decoration A'. This allows us to treat 
the overlap in an expansion that becomes asymptotically exact for large enough A^. 
In contrast to other large-N approaches|27l, this procedure occurs entirely within 
the class of SU(2) models. 

The scale of the spin gap is determined by the minimum of the {«,}. If this scale 
is sufficiently large, then perturbations of Eq.|4]will generate effective operators in 
the valence bond manifold. In particular, consider the perturbation: 

5// = 7 ^ S,- • Sj + ( (Si • Si J (S2 • Si, ) + (Si • S„j ) (S3 • s„3 ) j (5) 

where the first sum is over nearest neighbors and the second is over square plaque- 
ttes (the symmetric terms are not explicitly written). The degenerate perturbation 
theory involves accounting for the non-orthogonality of the valence bond manifold. 
In particular, if Sij = {i\j) is the overlap matrix, then we may consider the orthog- 

— 1/2 

onal basis \a) = Li'S'c,- |')- Because the overlap between different states is small, 
we may label | a) by its order unity component. In terms of this basis, the operator 
Eq.|5]becomes: 

Hap = {S-'l^5HS-'l^)ap (6) 

= l,{S-"^)ai{i\dH\j){S-'l'^)jp 
ij 

(7) 

where x = Uij is, a. matrix that is 1 if states |/) and \i) differ by the (minimal) 
length 4 loop and zero otherwise; n fi,i counts the number of flippable plaquettes in 
state lO; and r== -7x4(^+1). 

Therefore, up to small corrections, Eq.|5]acts like the RK quantum dimer Hamil- 
tonian in the orthogonalized basis and by decorating the lattice with a sufficient 
number of sites, the matrix elements beyond the dimer model can be made arbitrar- 
ily small. While this procedure will not capture the fine-tuned aspects of the dimer 
model, such as the v = f critical point for bipartite lattices, without taking ex- 
tremely large, the various gapped phases that appear can be obtained for a finite 
decoration. In particular, we may repeat this procedure on a decorated triangular 
lattice to obtain a model with a stable, SU(2)-invariant spin liquid phase. A similar 
construction may be used to realize the modulated phases discussed in the previous 
section. 
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We have shown that incommensurate structures and spin Hquids can, in principle, 
arise from purely local Hamiltonians that do not break any symmetries The useful- 
ness of quantum dimer models in answering this kind of question was hopefully 
conveyed. Along with the proofs of principle comes an understanding regarding 
physical mechanisms by which such structures may form. In the spin liquid ground 
states, the key elements were geometric frustration and strong quantum fluctuations 
through ring exchange. In the devil's staircase construction, the key ingredients were 
the fluctuating domain wall picture and competing interactions in a strong coupling 
limit. 

While the constructions involved rather complicated lattices and/or Hamiltoni- 
ans, it is likely that the ideas apply for simpler, though perhaps less (analytically) 
tractable, models where frustration enters in a less rigid way than a hard-core dimer 
constraint. For example, we may consider models with longer bonds ll28ll or doped 
dimer models ||29l . Therefore, we may hope the ideas presented here will provide a 
starting point for constructing physically realistic models with exotic phases. 
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